Absence of a boron isotope effect in the magnetic penetration depth of MgB2 
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The magnetic penetration depth A(0) in polycrystalline MgB2 for different boron isotopes 
^lOg^iig-j investigated by transverse field muon spin rotation. No boron isotope effect on 
the penetration depth A(0) was found within experimental error: AA(0)/A(0) = —0.8(8)%, suggest- 
ing that MgB2 is an adiabatic superconductor. This is in contrast to the substantial oxygen isotope 
effect on A(0) observed in cuprate high-temperature superconductors. 
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Since the discovery of superconductivity with transi- 
tion temperature Tc w 39 K in the binary intermetallic 
compound MgB2 , a large number of experimental and 
theoretical investigations were performed in order to ex- 
plain the mechanism and the origin of its remarkably high 
transition temperature. Experiments were done reveal- 
ing the important role played by the lattice excitations in 
this material m|. In particular, the substitution 

of the ^^B with ^°B has been demonstrated to shift Tc 
to higher temperatures 0, 0] j expected for a phonon 
mediated pairing mechanism. 

However, MgB2 differs from conventional superconduc- 
tors in several important aspects, including, for instance, 
the unusually high Tn and the anomalous specific heat 
. Calculations 0, H| based on the Eliashberg formal- 
ism support the experimental results |^ llOj , revealing 
MgB2 to be a two-band superconductor with two super- 
conducting gaps of different size, the larger one origi- 
nating from a 2D cr-band and the smaller one from a 
3D TT-band. The electronic cr-states are confined to the 
boron planes and couple very strongly to the in-plane vi- 
bration of the boron atoms {E2g phonon mode). This 
strong pairing, confined only to parts of the Fermi sur- 
face, is the principal contribution responsible for super- 
conductivity and mainly determines Tc- The 7r-statcs 
on the remaining parts of the Fermi surface form much 
weaker pairs. The double-gap structure explains most 
of the unusual physical properties of MgB2, such as the 
high critical temperature, the total Tc isotope-effect coef- 
ficient {a w 0.32 Qy, the temperature dependent specific 
heat , tunneling [3 and upper critical field anisotropy 

An interesting point to be clarified concerns the nature 
of the electron-lattice coupling. It was proposed 0, 0| 
that MgB^ is a non-adiabatic superconductor. Alexan- 
drov suggested that, because of the large couphng 
strength of the electrons to the E2£, phonon mode, the 
many-electron system is unstable and breaks down into a 



small polaron system, similar to the cuprate high temper- 
ature superconductors (HTSC) , where the charge carriers 
are trapped by local lattice distortions. Cappelluti et al. 
[isj proposed that the small value of the Fermi energy Ep 
of the a bands relative to the phonon energy ujph violates 
the adiabatic assumption [uph ^ Ep), opening up a non- 
adiabatic channel that enhances T,.. These models are in 
contrast to the conventional theory of superconductivity 
(Migdal adiabatic approximation), in which the density 
of states at the Ferixii level, the electron-phonon coupling 
constant, and the effective supercarrier mass to* are all 
independent of the mass M of the lattice atoms. Both 
non-adiabatic models explicitly predict in MgB2 a boron 
isotope effect (BIE) on the c harg e carrier effective mass 
TO* [3j ^3 ■ Similar models |l4i Il5| were already used 
to explain the large ox ygen isotope effe ct ( OIE) on m* 
observed in HTSC El 111 El [ll H [2l| by measur- 
ing the OIE on the magnetic field penetration depth A, 
a physical quantity directly related to the charge carrier 
effective mass. 

In this letter, a muon spin rotation (^SR) study of the 
magnetic penetration depth A(0) in polycrystalline MgB2 
for different boron isotopes (^''B/^^B) is reported. /j,SR is 
a powerful microscopic tool to measure the magnetic pen- 
etration depth A |22|. Indeed, in a polycrystalline type 
II superconductor with a perfect vortex lattice (VL) the 
average magnetic penetration depth A can be extracted 
from the muon-spin depolarization rate (j(T) oc A~^(T) 
[2^ . In our measurement, no BIE on A(0) was observed 
within experimental error [AA(0)/A(0) = -0.8(8)%], in 
contrast to the substantial OIE observed in cuprate 
HTSC El E3 El El El El. our results imply that 
polaronic or non-adiabatic effects in MgB2 are absent or 
negligibly small. 

The //SR experiments were performed on two polycrys- 
talline MgB2 samples containing ^-'^B (Mg^^B2) and ^°B 
(Mg^^B?). Full details of the sample synthesis are given 
in Refs [1 Ell ■ In brief, the two samples were synthe- 
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FIG. 1: Normalized field cooled (0.5 mT) magnetization as a 
function of temperature for Mg^''B2 and Mg^^B2 samples. 



sized using elemental Mg (99.9% pure in lump form) and 
isotopically pure boron (99.95% chemical purity, 99.5% 
isotope purity, <100 mesh) combined in a sealed Ta tube 
in a stoichiometric ratio. The Ta tube was then sealed 
in a quartz ampoule, placed in a 950°C box furnace for 
24 hours, and then removed and allowed to cool to room 
temperature. 

To examine the quality of the samples low field 
(0.5 mT, field-cooled) magnetization measurements were 
performed using a commercial Superconducting Quan- 
tum Interference Device. Figure ^ shows the temper- 
ature dependence of the magnetization for the Mg^^B2 
and Mg-'^°B2 samples in the vicinity of T^. The high 
quality of the two samples is revealed by the sharp tran- 
sition and the high extracted from the intercept of 
the linear extrapolations (Fig.P): T^i^^B) = 40.36(5) K, 
Tc("B) = 39.36(5) K. There is a clear isotope shift of 
ATc = Te("B)-Tc(i"B) = -1.00(7) K. The correspond- 
ing isotope effect coefficient as ~ —dln{Tc) / dln{MB) ~ 
0.29(2) (enrichment corrected) is in good agreement with 
previous results 0, ll| . 

The transverse-field ^SR experiments were performed 
at the Paul Scherrer Institute (PSI), Switzerland, using 
the 7rM3 /zSR facility. The samples used for the magne- 
tization measurements (see Fig. were pressed in disk- 
shaped pellets with 10 mm diameter and 3 mm thickness 
and cooled in an external magnetic field B^xt perpendic- 
ular to the muon spin polarization from well above 
to temperatures lower than T^. The measurements were 
taken in a field of B^xt = 0.6 T (the highest available 
at PSI), hi gh e nou gh t o avoid pinning induced distortion 
of the VL in miH. As shown in Fig. H for Mg"B2 
at two different temperatures, the local magnetic field 
distribution can be very well approximated by a single 
Gaussian, centered at a field lower than the external one. 
This again indicates the high quality of the samples and 
the absence of any normal conducting domains. From the 
width of the Gaussian field distribution, which is propor- 
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FIG. 2: Local magnetic field distribution, obtained from the 
Fourier transform of the muon spin precession signal, for 
Mg^^B2 at 2 (•) and 25 K (□). Solid lines are gaussian fits 
to the experimental data. The vertical solid line indicates the 
external field of 0.6 T. 

tional to the muon spin depolarization rate tj, the pene- 
tration depth A, that is the length scale of the variation 
of the magnetic field, can be extracted using the relation 

CX (T. 

In Fig. O the temperature dependence of a for the 
Mg"B2 (o) and Mgi°B2 (A) samples is shown. Below T^, 
a for both samples starts to increase and saturates at low 
temperatures T < 6K, in agreement with previous fiSK 
measurements jl^. The data for the two samples close 
to Tc show a clear isotope shift of ATc = — 1.2(2) K, in 
agreement with ATc deduced from the low field magneti- 
zation measurements (Fig. With decreasing tempera- 
ture, the values of a for the Mg-'^^B2 sample are systemat- 
ically lower than those for the Mg^''B2 sample. However, 
at low temperature they merge together, indicating that 
there is no substantial BIE on cr(0). 

In order to quantify this observation, we performed fits 
to the experimental data. It was suggested |25||, i26|| that 
for the two-gap superconductor MgB2, the temperature 
dependence of a can be written in the form: 

cr(r) = cr(0) - w ■ 5cr(Ai,r) - {1 ~ w) ■ da{A2,T) (1) 

with fc(A,r)=H^/„°° /(e,T).[l - fie,T)]de. 

Here, Ai and A2 are the zero temperature large and 
small gap, respectively, w is the relative contribution of 
the large gap to A~^(0), and f{£,T) is the Fermi distri- 
bution. For the temperature dependence of the gaps we 
used the conventional BCS A(T). In order to improve 
the ratio of the number of data point vs the number of 
fit parameters, the two gaps and w were considered as 
common fitting parameters for the two isotope data. As 
shown by the solid and dotted lines in Fig. the experi- 
mental data are well described by Eq. The fit yields: 
a{0f^ = 9.79(10)Ais-i, ^(0)'°^ = 9.95(ll)/^s-i, w ^ 
0.88(2), Ai = 4.9(1) and A2 = 1.1(3). AU these values 
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FIG. 3: Temperature dependence of a at Bcxt = 0.6 T for the two isotope samples Mg^''B2 (a) and Mg'^^B2 (o). The solid 
(Mg^^Ba) and dotted (Mg^Ba) lines are fits using Eq. (0. Inset: low-temperature region on a larger scale. The dotted and 
solid horizontal lines represent the weighed average values of (t(0) for T < 7.5 K for Mg^^Ba and Mg^^Ba, respectively. 



are in very good agreement with previous ^SR measure- 
ments performed by us on a natural boron MgB2 sample 
and by Ohishi et al. |2^. It is interesting to note that 
the high value of w implies that only a very small contri- 
bution to o-(O) originates from the 7r-band, in accordance 
with the experimental finding that the superfluid den- 
sity in the 7r-band is strongly suppressed by an external 
magnetic field 

The relative isotope shift of (t(0) is: (o'(O) ^ — 
a{0f^)/a{0f^ = A(7(0)/(t(0) = AA-2(o)/A-2(o) = 
— 1.6(1.5)%, corresponding to: 

AA(0)/A(0) -0.8(8)%. (2) 

For comparison, we calculated the relative isotope shift 
using a different and model independent procedure, tak- 
ing the weighed average of the experimental points for 
T < 7.5 K (see inset of Fig.OJl, where a{T) saturates. We 
obtained Acr(0)/CT(0) = AA-2(o)/A-2(0) = -0.5(8)%. 
Both results are compatible with zero BIE on the pene- 
tration depth A(0). 

Here, it is very important to recall that the two iso- 
tope samples used in the experiment were made with the 
same starting Mg for both the samples, and with ^°B 
and ^^B powders of the same mesh size (distribution of 
grain sizes) , and were synthesized under exactly the same 
conditions. Therefore, we can exclude any influence on a 
due to different grain size and to a difference in pinning 
or vortex dynamical effects. 

To check the reliability of our results, a second mea- 
surement on a set (set B) of samples from different source 
and preparation technique and with smaller Meissncr 
fraction, was performed in a field of 0.4 T. The results 
arc very similar to the first set (set A) shown above: 



AA~2(0)/A"2(o) = -1.5(1.7)% as compared to the above 
— 1.6(1.5)%. This shows that our result is intrinsic for 
MgB2 and holds for lower fields as well. A summary of 
the results for both sets of isotope samples is given in 
Table n Note that the values of cr(0)' ^ and cr(0)"'^ 
for set B measured in lower fields are larger than the 
corresponding values for set A. This is consistent with 
previous work |25l. E^. 

Theoretically, the zero temperature penetration depth 
is proportional to a density-of-states weighed average of 
a tensor involving the Fermi velocities. Detailed cal- 
culations within different formalisms have been carried 
out for MgB2 (see Refs. [13, H^)- For our purpose it is 
sufficient to use the simpler London approach consider- 
ing a free electron model and linking A(0) to the super- 
conducting charge carrier density and effective mass 
TO*, only considering different contributions from the tr 
and the tt bands. There is of course a direct connec- 
tion between the Fermi velocities and the effective mass 
(a band average), both of which are not bare quantities, 
but in general renormalized, e.g., due to coupling with 
the phonons. The London approach has the advantage 



TABLE I: Summary of the BIE results for a{0) obtained from 
the /iSR measurements of two sets of isotope samples. 









AA-2(o)/A-2(o) 




(/is-i) 


(/is-i) 




Set A 


9.95(11) 


9.79(10) 


-0.016(15)" -0.005(8)' 


Set B 


12.91(17) 


12.69(13) 


-0.015(17)" -0.016(30)' 



"from fit using Eq. 

''from low temperature average (inset Fig. UJ 
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of facilitating the comparison with theoretical predictions 
13, ll^ and results obtained on cuprate superconductors 
130 III III Hi 113, all of which are formulated within 
this approach. 

Unlike the cuprate superconductors with their ex- 
tremely short coherence lengths, MgB2 cannot be con- 
sidered as being in the superclean limit and we need to 
consider a possible impact of scattering. In a moderately 
clean superconductor the penetration depth is related to 
the effective mass m* by the following relation [2^ : 



l/\^ = [^,oe'|c\ns/m*)[l/{l + ^/l)], 



(3) 



where Us and m* are the superconducting charge carrier 
density and effective mass, respectively, ^ is the coher- 
ence length, and i is the mean free path. As already 
mentioned, the major contribution (~ 90%) to A^^ in our 
experimental conditions comes from the cr-band. There- 
fore Us, m*, ^, and £ in Eq.Q have to be considered 
as u-band values. It was estimated Isl I^ that in the 
(T-band {£,/£) a ~ 1/8, a value which is close to the clean 
limit (^/-^ <C 1). Therefore Eq. (jSJ may be approximated 
by l/A^ K [/xoeVc^]("-.s/m*). A shift in l/A^ due to the 
isotope substitution is then given by 



AA-2(0) _ An^ Am* 
A-2(0) ~ ~s TO^' 



(4) 



The contribution from the supercarrier density Ug is neg- 
ligible, as was alread y expe rimentally demonstrated in 
the case of HTSC d El El EH. Specifically, for MgBa, 
it can be argued that: i) by changing the isotope only the 
mass of the nuclei is changed and not the charge carrier 
density n. Furthermore, MgB2 is a stoichiometric com- 
pound; ii) x-ray diffraction measurements, performed on 
the samples used for the /^SR experiments, showed no 
substantial difference between the lattice parameters of 
Mg"B2 and Mgi"B2. This implies that the band struc- 
ture is not appreciably modified by the isotope substitu- 
tion. Therefore, assuming Aus/us ~ in Eq.Q) and ne- 
glecting the small 7r-band contribution, we can estimate 
the boron isotope effect on the a-band effective mass m* : 



Am* /m* 



-AA-2(0)/A-2(0) = 1.6(1.5)%. (5) 



Here we have used the value of the relative shift on A^^(O) 
obtained from the fit to Eq. There is no BIE on the 
(T-band effective mass within experimental error. 

Our result then suggests that non-adiabatic or pola- 
ronic effects in MgB2 are absent or negligibly small, and 
establishes an upper limit [Eq. ^ and Eg. lO] to any 
theoretical prediction of such effects 0, ^3 ■ This con- 
clusion is in contrast to cuprate superconductors, where 
a substantial oxygen isotope effect on m* was observed 
[lil ItH . nil 0, 12I I^ . Recent magnetization measure- 
ments on MgB2 under pressure |32l | show no pressure 



effect on the magnetic penetration depth A at low tem- 
perature, further supporting the main conclusion of the 
present work. 

In summary, /iSR experiments on polycrystalline 
Mg-'^''B2 and Mg^^B2 samples revealed no substantial 
boron isotope effect on the magnetic penetration depth 
at T = OK. From this finding we conclude that there is 
no substantial BIE on the effective mass m* of the charge 
carriers in the a band. This result suggests that MgB2 is 
a conventional phonon mediated superconductor without 
non-adiabatic or polaronic effects, in contrast to cuprate 
superconductors. 
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